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I . INTRODUCTION 



A. BACKGROUND AND RELATED WORK 

The study of electromagnetic energy transmission is but 
one important area in microwave and millimeter-wave engi- 
neering, where the electromagnetic waves are travelling 
through some transmission medium, which provides the link 
between the transmitting and receiving part of a transmis- 
sion system. 

In recent years, fin- line has gained in importance as a 
transmission medium in millimeter wave circuit constructions 
[Ref. 1-4]. Fin-line has been found superior to micro- 
strip at millimeter wavelengths as the former provides eased 
production tolerances, low dispersion, broad single mode 
bandwidth, moderate attenuation, better compatibility with 
hybrid devices, greater freedom from radiation and higher 
mode propagation, combined with the ability to construct 
simple transitions to conventional rectangular waveguide. 

Figure 1.1 shows a 3-dimensional view of fin- line. The 
stucture may be viewed as a slotline with a shield, a ridged 
waveguide with dielectric or a slab loaded waveguide with 
fins . 

The fin-line structure was first proposed for millimeter 
wave integrated circuits in 1974 by Meier [Ref. 1]. An 
early paper by Meier described the propagation mode as a 
variation of dominant mode in ridged waveguide. His proce- 
dure requires a test measurement to determine the equivalent 
dielectric constant of the fin-line structure. This is both 
expensive and time consuming. Knorr and Kuchler [Ref. 5] 
presented a frequency dependent hybrid-mode analysis of slot 
line with open boundary using the spectral domain technique 
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Figure 1.1 3 -Dimensional View of Fin-line Structure. 
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which was suggested by Itoh and Mittra [Ref. 6], in 1975. 
Subsequently, the paper of prime importance in the estab- 
lishment of spectral domain technique for analyzing the 
fin-line structure, so-called shielded slotline, was 
pesented by Knorr and Shayda in 1980 [Ref. 3]. One of the 
advantages of the spectral domain approach is that it is 
numerically more efficient than the conventional methods 
that work directly in the space domain. This is due prima- 
rily to the fact that the process of Fourier transformation 
of the coupled integral equations in the space domain yields 
a pair of algebraic equations in the spectral domain that 
are relatively easier to handle. Another important advantage 
is that the Green's function takes a much simpler form in 
the transform domain, as compared to the space domain where 
no convenient form of the Green's function is known to 
exist . 

B. PURPOSE 

In solving the open boundary slotline problem, only 
exponential functions arise and there are no numerical prob- 
lems during computation. For a closed boundary structure, 
however, hyperbolic functions are required. Knorr and Shayda 
[Ref. 3] found that overflow and underflow problems resulted 
during numerical computations. To eliminate these problems, 
extensive algebraic manipulation of the spectral domain 
equations is required. 

This thesis presents the method to calculate the charac- 
teristic impedance of fin- line without overflow and under- 
flow in equations. Therefore, this thesis is a direct 
extention of Knorr 's work. This work is accomplished by lots 
of algebraic manipulations. In particular, since hyperbolic 
sine and cosine functions cause overflow and underflow 
errors, the equations developed by Knorr need to be put in a 
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form where only the hyperbolic tangent function appears. The 
characteristic impedance is computed after the spectral 
domain technique to find the dispersion characteristic. 

Following the theoretical analysis, an explanation of 
the computer program used in determination of characteristic 
impedance is presented. Numerical results are then compared 
with known data for ridged waveguide, slab loaded waveguide, 
slot line and fin-line [Ref. 3]. This comparison estab- 
lishes the accuracy of the numerical results. 
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II. THEORETICAL ANALYSIS OF FIN LINE 



A. FIELD AND BOUNDARY CONDITIONS 

The fin-line supports a hybrid field. The axial compo- 
nants of TM and TE modes are then 



Ee = K? e r * (eqn 2.1) 

Hz = ICcftCXiW ^ (eqn 2.2) 



where the scalar potential functions **> ^ satisfy the 
Helmholz equation, and we assume lossless propagation so 
that f = £ jp . 

Further 




(eqn 2.3) 



2 . 2 . 

with kj, = , a .- | , 2, 3, defining spatial region of 

finline as stated in Figure 1.1. 

Through Maxwell's curl equations the transverse field 
components are then determined by these axial components and 
can be given as 



&=(rif ef2 

Hx=(r^ 
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(eqn 2.7) 




where propagation in the z direction is assumed. We will 
also assume here that £, = = £0 and £ 2 =&>£»* • 

At y = h i + D : 

Applying boundary conditions at the walls in region 1, 
tangential field components must be zero. 



At y = D : 

At the interface between region 1 and region 2, tangential 
field components must be continuous. 



Also the electric fields at y = D will exist only in the 
slot . 



Ea ta+P,.?) — o 



(eqn 2.8) 



Ex 1 0X,h, + Pi£) = ° 



(eqn 2.9) 



E 2 I (*,D,2) 



( eqn 2.10) 



txzC'X,P.2j) 



(eqn 2.11) 




(eqn 2.12) 




(eqn 2.13) 
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Similarly, tangential magnetic fields must be discontinuous 
by corresponding surface current densities. 






ix(>0 

O |Xl<x 



(eqn 2.14) 



Hxi (*,D,2 )-H)cl(>,c>/ 2) = j 



> VZ 

0 ? W e 
o 



IXUf 

l*Kx 



(eqn 2.15) 



AT y = 0 : 

The tangential field components at the interface between 
region 2 and 3 must also be continuous. 

E22 0 *,o,*) = Ez 3 C^)° 1 ^ 3 

E xz(*jO>£)- Ex3(oc ,o t z) 

Wl2.(.*,0,£)=- C +>0,2) 

Hx2_(*,o ,2) ^ 

At y = -h^ : 

Once again at the shield wall in region 3, the tangential 
field components must be zero. 



E£3C'X,“hz,£-) =0 


(eqn 2.20) 


£^3 COC , "lii; B) — O 


(eqn 2.21) 



(eqn 2.16) 
(eqn 2.17) 
(eqn 2.18) 
(eqn 2.19) 
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At x = t.b/2 : 

The final boundary conditions occur at x = b/2 where the 

tangential components must be zero in all regions. 



E*; 

=0 



( eqn 2.22) 
( eqn 2.23) 



B. SPECTRAL DOMAIN APPROACH TO DISPERSION CHARACTERISTC 

The scalar potential functions can be transformed into 
the spectral domain via Fourier transform defined as 



h 



= §GM) = j” ^ M ■ 



(eqn 2.24) 



The scalar potential functions satisfy Helmholz’s equations 
in the three spatial regions, thus 



Vt fa t I Cci ~ ° 



(eqn 2.25) 



where denotes the two dimensional Laplacian operator in 
the transverse (x,y) direction. The Helmholtz equation 2.25 
are transformed into 






3 



(eqn 2.26) 



where it- £ - + p - K . Above equation has solu- 
tions after applying boundary conditions at y = D + hj , and 
y = -h 2 . 






(eqn 2.27) 
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£2 6m) - Vij + c%^)cosU 


(eqn 2.28) 


§3 (*14) = D'CaVMrvMk ) 


(eqn 2.29) 


= frWcesU )>, (D+l.,-3) 


(eqn 2.30) 


f z <•&«, J >ijj + c l (°U-)cosU 


(eqn 2.31) 


l3(^-3) = D h (o^)c»sh 


(eqn 2.32) 



where 




rtlTT 

b 

b 



^ ei/e/rt 




( eqn 2.33) 



(eqn 2.34) 



The following observation about the solutions in region 
2 must be emphasized. Any wave on this inhomogeneous wave- 
guide structure will partly travel through air and partly 
through the dielectric slab. It is important to observe at 
this point that may be less than zero in any of the 
three regions of the structure under certain conditions. 
When o6n = 0 and approaches k 0 (where k<> is the wave 
number for free space), is less than k^ and so Vi < 0 • 
Under this condition the hyperbolic functions in all three 
regions are replaced by trigonometric functions. If k© < $ < 
k 2 , then ^ and are greater than zero and < 0 for some 
values of and the trigonometric functions replace the 
hyperbolic functions in the spatial region 2 only. This 
suggests that the nature of the field is dependent upon the 
values of the transform variable, . For the conditions 
when < 0, y*| replaces such that (V/. ) = - yT* The 
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eight unknown 
dent , but can 
nuity conditions 
between the three 
If we denote 
current density 

&&0 = F { e*(oj 

BtCU) = F [ ejW)} 

* F [ j,(-X)j 

3z«0 = F {>(*)} . 

The resulting set of linear 
matrix form as follows: 



are not indepen- 
through the conti- 
at the interfaces 

z-directed 

( eqn 2.35) 
(eqn 2.36) 
(eqn 2.37) 
(eqn 2.38) 
equations may be written in 



coefficients A through D 
be related to each other 
of the field components 
spatial regions. 

the Fourier transforms of x- and 
and electric field component by 





A e ' 

B* 

• 

• 


. 


i — 

O 0 • • 




* 

k u > 




« 

Ex 

. S - 








/ 

0 o 
\ 


K) 


• 

• 

* 

C* 

D h 

k > 




• 

\ 

£ 



(eqn 2.40) 
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The matrices ([ME] and (MJ) differ in only the last two rows. 
Each is a square 8x8 matrix. Using equations (2.39) and 
(2.40), we may write 






' O ' 




' o ' 


o 




o 


• 




• 


# 




• 






r>~s 


Ex 






/V 






Ez 




3k 


k J 




^ ) 



(eqn 2.41) 



From equation (2.41), using the four elements in the lower 
right hand corner of the matrix MJ ME , we obtain 





tf2_ 




' &(*.)' 




'J»w' 


* 






> 




JzM 



where the 2 x 2 matrix [g] contains the Fourier transforms 
of the components of the dyadic Green's function for this 
structure . 

A solution to equation (2.42) is obtained using the 
Method of Moments [Ref. 7]. For this problem, we have 
chosen to approximate the field between the fins as shown in 
Figure 2.1: 



e*(*) = 



f i 
. 0 



1*1 i X 

elseiOkere 



(eqn 2.43) 



e*(0 = o 



(eqn 2.44) 
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Figure 2.1 Assumed Electric Field Component in Slot 
in x-direction versus x for Fin-line. 
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The dispersion problem is now reduced to the form 






(eqn 2.45) 



where 

%440-jt 



SC/A (o^n h)^2_) 

oCr\ 'd/x 



A numerical search for the value of which satisfies equa- 
tion (2.45) yields the propagation constant for the dominant 
fin- line mode. 

From the equation (2.39) 



/Vilii f\ e + /Witz B e 4- /*Wi 3 C e = O 


( eqn 2.46) 


AHii pf ■+■ 4 Ma 3 *4“ 




A^i5 A* 1 4 /v*U& B*' 4 /VM d = 0 


(eqn 2.47) 


wfl A e =• 


(eqn 2.48) 


4 A* 1 = p 1- Ex 


(eqn 2.49) 


*133 <C e -V- /ma4P e - - 0 


(eqn 2.50) 


m 43 c e 4 /m 44 if + mu B 1 ' + M 48 P 1 ' - 0 


(eqn 2.51) 


/i w 5TjC t ' + /WVcrS [f = 0 


(eqn 2.52) 


B € 4- /WU4.D e C* 1 d- Mi 8 D* 1 = O 


(eqn 2.53) 



From equations (2.50) and (2.52) 
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(eqn 2.54) 



/*33 c e- 




( eqn 2.55) 



Using equations (2.48) and (2.49) 




(eqn 2.56) 




( eqn 2.57) 



Substituting in equations (2.51) and (2.53), we obtain 



* ( eqn 2.58) 

-t, _ //W44/*k3 _ SW ) c e + {M4*^ \ c ly 

* W^/^34 U*t4£/*SS/ * (eqn 2.59) 

Substituting equations (2.48), (2.49), (2.51), d (2.53) in 
equations (2.46) and (2.47) we obtain 




( eqn 2.60) 




(eqn 2.61) 
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where 



[D[ — 4|l - Oil d.12. 

— j (4(1 chi — y±p) • 



We define 



a„= U 8 +Mta(^-fiL)l 

*- \/VHX /M34/J 

= otn O^p} (eqn 2.62) 

*• /BUi -> 

- j d|2. (V^P) S|7 mK CVaP) (eqn 2.63) 



r- fMlA/fi 133^V . 

«« = Ki+^-u fe/WH j + 



/ywi44-/^33 


y^43 -\-j 


\/*U4 /1M34. 


/vM4t yJ 



= dUl C^lP) 



(eqn 2.64) 



(^1 - [^/IMay *t* 




i-y*Vbx 



f^Sj^L 
V/M61 /W5-8 




= J o(xL O^P) SwH ^ C^-P3 



(eqn 2.65) 



Also we define 
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dii = -(Kopfp + 



OJ£iD)CLv) (fejPf 
CtiP) OLdJ 



C^jbwihCC&PX^}] 



] 






(eqn 2 



( eqn 2 



4,_ d-uoXfP) L cutiP) (tp) CkiiP) 3 - 

(W^jlCM (fccrpf -, 
Cy^AvwKCO^X 11 ^ + 



( eqn 2 



4* = ^* m * m ^* 

L GfoD)XtoirvrthOlPj) 

r 6<MDfC^>r L - i rM _ i'll 
^ iGdmcu&Ktoyvl (biof- }) . 



(eqn 2 



Nomalized constants are defined as follows 



Cfopf = («fCi -C%.T) (p/a) 2 - 

OiiP)^ =• C iT f- C^r - CtytfYJ ( p /a/- 

U)MO = -2-40 IT 2 - (t>/70 

^£,D = ^ C°M) 



. 66 ) 



. 67) 



. 68 ) 



.69) 
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i-yji °OVa) 

= iir C p A)CVX) 

/« 2IT (jb ) 4> h eitert 

. OH)TT (-§-) <t )l ' 

CV,Df= CV»of = C^>of+ (arfCCWf-O C D Af 
<&*■■= (<4,d/+ (aTfCcV^f-^J C D Af . 



C. CHARACTERISTIC IMPEDANCE 



The definition of the characteristic impedance for an 
ideal TEM transmission line is uniquely given by static 
quantities. Since the fin- line supports a hybrid mode, no 
unique definition of the characteristic impedance can be 
found . 

A useful definition, however, is 



£> = 



2Pc 



V 1 



where V 0 is the slot voltage defined as 



Vo =• 0^ ^ ^ = 1 • 



(eqn 2.70) 



( eqn 2.71) 
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1 . 



e^(x) is arbitrarily selected as 1/W so that W'’'e*(x) = 

P av <j is given by 



P*v/^ — \ Re E 0 & 

nS 



( eqn 2.72) 



Pasval’s theorem is applied to eq(2.88) to obtain 

p«3 = ± & i £ * CEx 64., a ) Hf (56., -j ) 



This expression must be evaluated in each of the three 
regions of the fin- line shown in Figure 1.1 Therefore the 
power flow may be expressed by 



H = i ir- t / pl + P 2 + 



(eqn 2.73) 



Since P«\/| can be determined after finding the value of 
. Equations including only hyperbolic tangent functions 
and the slot voltage V 0 have been developed the character- 
istic impedance. The lengthy algebraic manipulations are 
shown in appendix B. 
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III. COMPUTER PROGRAMMING 



A. NUMERICAL ANALYSIS 

The computation of characteristic impedance is based 
upon the solution to the dispersion chracteristic problem 
for the fin- line under consideration. In other words, the 
wave propagation constant or wavelength ratio X. /?\ must be 
known before any other investigations can be started since 
only in this case are the scalar potential functions in the 
transform domain known. The computer program for the wavel- 
ength ratio X I ‘K is already developed by Prof. Knorr. The 
next task is the preparation of the appropriate equations 
for the numerical evaluation of the time average power flow. 
For ease in numerical calculations and for programming 
purposes all geometric parameters are normalized as follows; 
hj/D ; fin location relative to the positive "y" side 
wall normalized with respect to D 
h 2 /D ; fin location relative to the negative "y" side 
wall normalized with respect to D 
b/D ; waveguide height normalized with respect to D. 

It is observed that for the power flow computations the 
infinite integration in equation (2.88) is replaced by an 
infinite series due to the discrete nature of the transform 
variable JLm . Preliminary numerical investigations of the 
coefficients of these series indicated an even distribution 
with respect to the variable so that computation time can 
be saved by summing over a half interval only. 

The time average power flow equations are prepared for 
six different cases as follows; 

case 1 ; (V ( D f~ < 0 in region 1 

case 2 ; (V, D) 2 " > 0 in region 1 

case 3 ; (V^Df* < 0 in region 2 and (Y3D f~ < 0 
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(V,D f > 
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case 
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(&D f 
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in 
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2 
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'S 


0 
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case 


6 




O^Df 


> 


0 


in 


region 


2 


and 


(V»D f > 


0 






in 


region 


3 

















After this preparation of the equations for the computer 
programming, the only remaining task, is to investigate the 
numerical . integration with regard to its limits and the 
procedure to be used. One of the limits of the summation 
discussed previously extends to infinity so it is necessary 
to determine an appropriate at which to truncate the 
computations . 

From the representative examples shown in Figure 3.1, it 
is seen that the coefficients in this infinite series are 
found to decay rapidly so that a finite approximation yields 
good results. Figure 3.1 shows the characteristic impedance 
as a function of n, the number of terms in the truncated 
series, for various slot widths. The value n = 50 is 
sufficiently large to be considered infinite for all prac- 
tical purposes. It is noted that when w/b = 1, Ex(j^) = 0 
for n > 0. Therefore, whenever w/b = 1 only the n= 0 term is 
computed in impedance calculations. 

B. COMPUTER PROGRAMMING AND LIMITATION 

The first part of the FINIMP computer program finds the 

/ 

A/A which makes equation (2.45) equal to zero. The assump- 
tion is that A/A is between 0.1 and 3 in the giga hertz 
range and continuously decreases as the frequecy is 
increased. Since this program is very sensitive there is a 
possibility that the wrong value is searched. In that case, 
we can easily notice because there is abruptly jumped value. 
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CONVERGENCE TEST 



' 



o 




Figure 3.1 Characteristic Impedance Z 
for a Fin-line with b/D=18.8 h|/D=18.8 
D= .005 £>-=2.2 f = 40 . OGHZ . . 



vs . Iteration 
hz/D=17.8 
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At that time, if the assumed range of X/A i s resticted to a 
smaller interval than before the correct value can be 
obtained. 

The rest of the FINIMP computer program is the computa- 
tion of characteristic impedance. Since all equations are 
composed of hyperbolic tangent functions, there is no over- 
flow in the computer program. But if (V/.D) 2 " < 0 tangent 
functions replace the hyperbolic tangent functions. At that 
time, there is a possibity of overflow because tan n IT / 2 is 

infinite for n= 1, 3, 5, In that case the tangent 

function should be set to some value which is the maximum 
value of computer ability (46 6? - (1-16 - ^) > | tan(x) | ) . 
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IV. NUMERICAL RESULTS AND COMPARISONS 



To check the accuracy of the numerical results gererated 
by the computer program, comparisons are made with data 
available in the literature for ridged waveguide, slab 
loaded waveguide, slotline, empty waveguide and fin- line as 
outlined in [Ref. 3]. 

A. RIDGED WAVEGUIDE 

When = 1 and w/D < 1 or when the dielectric substrate 
thickness D is reduced to zero, the fin- line structure 
becomes ridged waveguide with zero thickness ridges. 

The impedance of the ridged waveguide has the following 
relations [Ref. 8], 



z 



o 



2?o fco 
0 - 



( eqn 4.1) 



When the guide width and the slot width are equal, the 
ridged waveguide becomes an ordinary rectangular waveguide 
for which 
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2A_ 

b 



(eqn 4.2) 



When the guide width is increased at fixed slot width, the 
characteristic impedance at infinite frequency and the free 
space cutoff wavelength go asymptotically to 
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Figure 4.1 Characteristic Impedance Z vs. 
Frequency for a Ridged Waveguide. 
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Figure 4.2 Characteristic Impedance Z vs. 
Frequency for a Ridged Waveguide. 
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v Im* .ppL 



V trU a 



( eqn 4.3) 



The cutoff wavelength ^ was determined from [Ref. 11] for 
several values of W/b. 

Equation (4.1) was then used to calculate the character- 
istic impedances and these values were compared with the 
results using spectral domain method of the computer program 
discussed in chapter JL . Two results are very close 
together. Figure 4.1 and 4.2 shows the computer results. 
These figures agree with the previous Knorr and Shayda's 
results . 

B. DIELECTRIC SLAB LOADED WAVEGUIDE 

A variation of the fin- line structure of Figure 1 where 
w/b = 1 and Cy > 1 results in a dielectric slab loaded 
rectangular waveguide. The slab loaded waveguide has also 
been studied by Vartanian et al [Ref. 9]. They consider a 
guide with the slab centered in the waveguide and they 
present an analytical expression for the voltage impedance 
Zpv at the center of the slab. The impedance computed in 
this analysis was specified at the edge of the dielectric 
slab. Therefore a relationship between the two impedances 
must be defined before any comparison can be made. The impe- 
dance at the edge of the slab can be easily defined by the 
expression 




(eqn 4.4) 
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SLAB LOADED WAVEGUIDE 




Figure 4.3 Characteristic Impedance Z vs. 
Frequency for a Slab Loaded Waveguide. 
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Figure 4.4 Characteristic Impedance Z vs. 
Frequency for a Slot Line Waveguide. 
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where is the field at the edge of the slab and 

is the field at the center of slab, thus 



- Zpv c<?s2 "GH) 



(eqn 4.5) 



where the various quantities are defined in [Ref. 9] as 



5 = 




(eqn 4.6) 



(40 _ ^ G.) “(IjO (eqn 4.7) 

— ( 23 9 i ('^‘) r^-C^J (eqn 4.8) 

and c = D is the dielectric slab thickness. 

The characteristic impedance of a slab loaded guide is 
computed using the spectral domain method. The results were 
compared with equation (4.5) results and [Ref. 3]. Good 
agreements are obtained. Figure 4.3 shows the slab loaded 
wave guide impedance . 

C. SLOTLINE 

If w/D < 2 and £y" is sufficiently high for the fin- line 
structure of Figure 1.1, the field is tightly bound to the 
slot. For this condition the presence of the shield will 
have little effect if the walls are sufficiently far removed 
from the slot. In this case the fin-line will behave like a 
slotline. This behavior is illustrated in Figure 4.4 where 
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the characteristic impedance of a fin- line with w/D =1, £ r 
= 20 have been plotted. Also these results are in good 
agreements with [Ref. 3]. 

D. SEVERAL FIN LINE IMPEDANCE CURVES 

Fin-lines are generally enclosed with a shield that is 
compatible with the dimensions of the standard rectangular 
waveguides for the millimeter wavebands. Above 22GHZ , all 
these guides have aspect ratios b/a = 0.5. Further, the 
fins are most often centered in the guide and printed using 
D = 0.005 inch substrate with £y = 2.2. 

Figures 4.5 - 4.12 show the impedance computed for fin- 
line with WR(28), WR(19) and WR(12) rectangular waveguide 
shields. These figures may be compared with data for the 
same structures as presented in [Ref. 3]. Such a comparison 
shows excellent agreement for all but the small values of 
W/b. For small values of W/b (W/b = 0.1, W/b = 0.05, W/b = 
0.02), the results of the FINIMP program show lower impe- 
dances than the results presented in [Ref. 3]. As the W/b 
is smaller, the difference is larger; for W/b =0.1, 2-5 
ohm is lower, for W/b =0.05, 9 - 15 ohm is lower and for 
W/b = 0.02, 19 - 25 ohm is lower. 
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FIN-LINE WITH WR(28) SHIELD 



r 



© 




Figure 4.5 Characteristic Impedance Z vs. Frequency 
for a Fin-line With b/D=28.0 h|/D = 28.0 h 2 ,/D = 27.0 

D= .005" y. = 2 . 2 . / 
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FIN-LINE WITH WR(28) SHIELD 




Figure 4.6 Characteris c Impedance Z vs. Frequency 
for a Fin-line With b J=28.0 h,/D=28.0 h z /D=27.0 

D= . 005 £ r =2.2. 
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IMPEDANCE(OHM) 
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FIN-LINE WITH WR(19) SHIELD 




FREQUENCY (GHZ) 



Figure 4.7 Characteristic Impedance Z vs. Frequency 
for a Fin-line With b/D=18.8 h|/D=18.8 h 2 /D=l /.8 

D = . 005" ty =2.2 . 
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FIN-LINE WITH WR(19) SHIELD 




LEGEND 
o W/B=0.1 
o W/B=0.( 
W/B=0.02 



46 48 50 52 54 

FREQUENCY (GHZ) 



Figure 4.8 Characteristic Impedance Z vs. Frequency 
for a Fin-line With b/D=18.o hj/D=18.8 ha/D=l/.8 

D= .005” ty =2.2. 
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FIN-LINE WITH WR(12) SHIELD 
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Figure 4.9 Characteristic Impedance Z vs. Frequency 
for a Fin-line With b/D=12.2 h|/D=12.2 h*/D=lI.2 

D= .005" lr =2.21 
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FIN-LINE WITH WR(12) SHIELD 




Figure 4.10 Characteristic Impedance Z vs. Frequency 
for a Fin-line With b/D=12.2 hi/D=12.2 hz/D=lI.2 

D= . 00d £ r = 2.2. 
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FIN-LINE WITH WR(19) SHIELD 



o 




Figure 4.11 Characteristic Impedance Z 
for a Fin-line With b/D=18.8 h|/D=28.2 

D= . 005" £ r =2.2. 



vs . F requency 
hz/D=8 .4 
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FIN-LINE WITH WR(19) SHIELD 
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Figure 4.12 Characteristic Impedance Z vs. Frequency 
for a Fin-line With b/D=18.8 h|/D=28.2 hi/D=8.4 

D= .005" £ y =2.2. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The spectral domain technique used in conjunction with 
Galerkin’s method has been presented to calculate the char- 
acteristic impedance for the dominant mode of the fin- line. 
It has been shown that a matrix formulation of the problem 
permits the elements of the dyadic Green's function to be 
calculated. Solving the matrix equations leads to 
containing only hyperbolic tangent function. These equa- 
tions circumvent the overflow and underflow problems which 
occurred using the previous formulation presented by Knorr 
and Shayda. 

Numerical results obtained using this method have been 
presented and compared to other existing data. Good agree- 
ment has been obtained in all cases thus establishing the 
accuracy and applicability of the method for the full range 
of structure parameters. 

There is a possibility that tangent functions may cause 
an overflow problem if (VaO)*^ is less than zero. In this 
case, the value of tangent function can be obtained within 
the capability of the available IBM 3033 used. 

In this thesis, paticular interest is devoted to the 
computation of fin-line impedance. Fin-line may exhibit the 
characteristics of ridged waveguide, dielectric slab loaded 
waveguide, slot lines, and conventional rectangular wave- 
guide. All of these structures are fin-line substructures. 
So, the computation of fin- line impedance permits all of 
these structures to be analyzed. 
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B. RECOMMENDATIONS 



Coupled fin-line will find future use in building direc- 
tional couplers and filters. For this purpose, the normal 
mode wavelengths and impedances are required. A special but 
important case is that of symmetrical lines. Coupled fin- 
lines for which the normal modes are odd and even. The 
program described here may be extended to cover this case by 
following the procedures outlined by Knorr and Kuchler for 
coupled slotlines [Ref. 5]. This should be accomplished to 
improve the utility of the program described in this thesis. 
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APPENDIX A 



SPECTRAL DOMAIN MATRICES 

The continuity conditions are transformed via equation 
(2.24) into the two dimesional Fourier domain. The solutions 
to the two Helmholtz equations given by equations 
( 2 . 27 ) - (2 . 32 ) are substituted. Finally a matrix form of 
linear equations is derived as follows; 
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The matrix elements of fMg} and [Mjl are normalized at 
this point with respect to D, the dielectric substrate 
thickness. The normalized matrix elements are now presented 
in two forms. The first of the element equations is for 
(Ya D) 2 * > 0 and the second is for (l/i D) 2 ' < 0. For the 

matrices [Mel and CMj 3 , the elements m ( | through m 68 are the 
same . 
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APPENDIX B 



TIME AVERAGE POWER FLOW 



The following expression for the coefficients A 
are derived from equations (2.54) - (2.61). 

If (V,D) Z < 0 
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If ( VxD f~ > 0 and (V 3 Df< 0 
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REGION 1 

The term Pj^ will be used for the power flow in region 1 
for the case (y ( D) 2- > 0. 

P,a = ' i(p)Ke £ „ { + Cp>)<W> 
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For the case ( V, D) 2- <0, Vi is imaginary in which case 
OV'D) 2 - = -(V|D) X . The power flow, , for this case is 
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CfP)(UAD) 
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+ (p>f(s ^ 



04 , 0 ) ] [ 2 Ct'o)C%) + sr* [iCYi'-t-Xto]]] 



REGION 2 

For region 2, the power flow expressions are the same as 
for unshielded slotline except that the Fourier integral is 
replaced by a summation and the interval 2 is replaced by 
the interval b. 

For the case (VzD) 2, >0, Vt is real and 
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For the case 
OzDf = -(Va-Df. 



(Va-D f < 0, Vi is imaginary in which case 
The power flow, , for this case is 
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REGION 3 

For the case (V 3 d f > 0, is real and the power flow, 

P3<\> is 
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For the case (V3D) 2 " < 0 , V3 is imaginary in which case 
(V^D) 2- = -(VjD) 3 *. The power flow, j for this case is 
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APPENDIX C 



COMPUTER PROGRAM 'FINIMP' 
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